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ABSTRACT 



This paper presents a new numerical method for solving 
Schroedinger ! s equation in radial form in order to get a set 
of wave functions for use as starting values in the solution 
of the Fock equation,, As a necessary adjunct to the solution, 
an improved numerical method is given for calculating potentials 
for use in the Schroedinger equation. Finally, a method is 
outlined for solving the Fock equation numerically, using the 
wave functions obtained above as starting values. 
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ENERGY OF THE MOST LOOSELY BOUT® ELECTRON 



FOR IONS OF COPPER, SILVER, AND GOLD 



If the energy and distance from the nucleus of the most 
loosely bound electron in an ion of a metallic lattice can be 
found, it is possible to use this data to calculate certain 
physical constants of the metal, e.g. the elastic constants. 

Wigner and Seitz ^ have shown how to find this binding 
energy, neglecting exchange forces, by solving Schroedinger * s 
equation in radial form: 

(D -ii c£5 [vW- Oe - o 

N 1 2rw\ A. r 



for various values of E. It is found that for each such value 
of E there are several radii at which the condition 

iiL - JL 

dr r 

is satisfied. Most of these points are very close to the 
nucleus and are essentially the same for all values of E. If 
the outer two radii are plotted against E a curve results 
which exhibits a minimum. It is this minimum which gives the 
required data. 

p 

Fuchs has extended the procedure to take into account 

3 

exchange forces. His method involves solving the Foci equation: 



^E. Wigner and Seitz, F. , Phys. Rev., 4^, 804 (1953) 
2 K, Fuchs, Proc. Roy. Soc., A, 1^1, 585 (1955) 
h, Fock, Zeits. F. Physik, 81, 195 (1955) 
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( 5 ) 



-%1 d^B +[v(<)-e]E- 0 . 

2. crs\ 0 

The most practical method for solving this equation is 
first to solve equation (l) and then use the wave functions 
obtained to evaluate the integral in (5)* This makes the right 
side a function of the independent variable instead of zero, 
When this resulting equation is solved the new wave functions 
are then used to re-evaluate the integral and the procedure is 
repeated as required* 

As the first approximation, then, a solution of (1) is 
needed* To get it, it is necessary first to calculate the 
values of V to be used* 

From the work of Hartree and others , tables of values of 
Z, the effective charge, are available as functions of (O , the 
radius measured from the surface of the nucleus. These values 
have been accepted in this paper, although there may be errors 
in them since most were calculated by hand before the era of 
the electronic computer. The required values of v, however, 
were recalculated, where 




and Z is the effective nuclear charge for potential which is 

p 

different fron Z. The values of v were calculated by 

^D. Hartree, Proc. Roy. Soc., A, l4l , 282 (1955) 

A. Douglas, Hartree D., and Runciman, W. Proc. Cambridge 
Phil. Soc., £1, 486 (1955) 

M. Black, Manchester Memoirs, 79 . 29 (195^0 
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integrating the field intensity 



(5) ' ) - V* 



n 



The usual method^ of calculating v is to solve the 
differential equation 

d-^£ _ ’ip - ^ 



(6) 



**■ e 



e 



by Suler*s method, which is simple but inaccurate. It can be 
shown 0 that the error is of the order of 



(7) 



s 



(f.-fO 1 " (A 






Some simple algebra and substitution reduces this expression to 

(s) e = (JdQ \i(e.) -~-(e,)\ 

zp, 

since Z is assumed to be linear between any two adjacent radii. 
Graph 1 gives the relative % error for Cu + resulting from this 
calculation. It is large enough so that a more sophisticated 
method of solution, equation (5), was decided upon. The 
evaluation was done using a Gaussian p“P oint quadrature formula, 
which is based on the assumption that between any two neighboring 
radii Z is approximated by a third degree polynomial. 



To evaluate the integral in (5) > 




5d. Hartree, Proc. Cambridge Phil. Soc., 24, ll4 (1928) 

^Kunz, K., Numerical Analysis, McGraw-Hill, N.Y., 1957* pgs 177-8 
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it is first necessary to normalize the limits to the range 0 to 1 
in order to use readily available values of the Gaussian coef- 
ficients . This can be done by the substitution 



( 10 ) 



u. - 1_ 



x r - (f>, +-f_) 



Equation (9) then becomes 

, , / (a (f v -p,) + (, 

(11) A ({0- (2H.-,)({(.-p 1 )+ (Pi+'P.)’ - 






The coefficients for the Gauss 5-point formula are 
it 0 * . II2.-7 016,7 
= • 5 

Thus, (11) becomes 



(12) 



, S_- 

Vvo r 

w,-- * 
£ 




•+ 









r< =. \ . i-74s<3\664 ^ (2>- <1-540 2=34 

It can be shovrn^ that the error is of the order of 

05) £ -- b ^ W-nA(0 

L u 0 

where f represents the integrand in (9)« The magnitude of the 
relative error for Cu + is plotted in Graph 1, It is never 
greater than . 25 % 

To use (15) it is necessary to develop an interpolation 

^Lanczos, C., Applied Analysis, Prentice Hall, N.J., 1956, 

pg3 bob-7 
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formula for calculating the 5 numerators on the right from the 
available table of values of Z* This has been done using a 
Gregory-Newton second degree divided difference interpolation 
formula,, A third degree formula which would be consonant with 
the basic assumption involved in the Gaussian method 

would have increased the labor of programming significantly 
and would have made no difference in the final result, to the 
limits of accuracy retained* 

The resulting threo formulas arfi of tho form 



Appendix 1 gives the CDC~l6o4 program for calculating v, 



- . 18181818181818. The first v calculated by the program is 
that for ^ =5 Bohr radii. The values obtained are given in 
Table 1 . 

Equation (l) can now be solved. It is first desireable to 
write it in terms of Bohr radii and Rydberg energy units to get 
a dimensionless equation. 

Let 





calculated by equation 




(17) 




where a is a Bohr radius. Then 



o 



(18) 



i dd-e 



d. <Xi olp 
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TABLE I 



Cu + Ag + Au + 



r 


✓ — 
Z 




“ *> 

ur 


S' 

Z 




XT 


/- — " 

z 


-u" 


.000 


29.000 






47.00 






79.00 




.002 








47.00 


- 


15910. 






-3" 

o 

o 

• 








46.99 


- 


7458. 






.005 


28.992 


- 


5668.2 












.006 








46.95 


- 


4987. 






.008 








46.91 


- 


5705. 






.010 


28.948 


- 


2769.7 


46.85 


- 


2919. 






.015 


28.875 


- 


1805.0 


45.65 


- 


1945. 






.020 


28.759 


- 


1524.4 


46.56 


- 


1458. 


77.55 


- 5502. 


.025 


28.615 


- 


1057.5 


46.o6 


- 


1129. 






.050 


28.451 


- 


847.07 


45.78 


- 


920.6 






C55 


28.280 


- 


711.96 


45.50 


- 


771.6 






.04 


28.106 


- 


611.27 


45.25 


- 


659.8 


75.55 


- 1575. 


.05 


27.765 


- 


471.60 


44.75 


- 


509.6 






.06 


27.451 


- 


579.58 


44.20 


- 


409.2 


72.86 


- 955.5 


.07 


27.161 


- 


514.57 


45.62 


- 


557.8 






.08 


26.884 


- 


266.50 


42.5*5 


- 


284.7 


70. i4 


- 655.6 


.09 


26.608 


- 


229.14 


42.26 


- 


245.8 






.10 


26.520 


- 


199.71 


41.52 


- 


211.6 


67.76 


- 485.1 


.11 








40.78 


~ 


185.7 






.12 


25.689 


- 


156.27 


40.04 


- 


164.5 


65.41 


- 571.9 


.14 


24.977 


- 


126.05 


58.71 


- 


152.8 


62.92 


- 295.5 


.16 


24.209 


- 


104.07 


57.55 


- 


109.5 


60.51 


- 24o.5 


.18 


25.415 


- 


87.555 


56.48 


- 


91.81 


57.75 


- 199.4 



7 



TABLE I (Continued) 





Cu + 






A S + 






Au 


+ 




r 


' Z " 




IT 


•"z 




IT 


'IT 




N 

IT 


.20 


22.646 


- 


74.745 


55.48 


- 


11.99 


55.44 


- 


168.0 


.22 


21 .929 


- 


64.625 


54.48 


- 


66.91 


55.40 


- 


145.5 


.24 


21.279 


- 


56.452 


55.45 


- 


58.07 


51.66 


- 


125.4 


.26 


20.699 


- 


49.755 


52.51 


- 


50.78 


50.09 


- 


107.1 


.28 


20.187 


- 


44.124 


51.16 


- 


44.76 


48.62 


- 


95.55 


.50 


19.756 


- 


59.579 


29.98 


- 


59.75 


47.15 


- 


82.14 


.52 








28.80 


- 


55.55 


45.62 


- 


12 . hi 


.54 








27.66 


- 


51.95 


44.05 




6 4.21 


.55 


18.791 


- 


50.224 














• 56 








26.57 


- 


28.89 


42.56 


- 


57.15 


.58 








25.55 


- 


26.24 


4o.65 


- 


51 08 


.4o 


17.964 


- 


25.650 


24. 61 


- 


25.96 


58.94 


- 


45.84 


.45 


17.114 


- 


18.765 


22.65 


- 


19.48 


54.75 


- 


55.65 


.50 


16.195 


- 


15.050 


21 .16 


- 


16.06 


50.96 


- 


28.41 


.55 


15.196 


- 


12.189 


20. 04 


- 


15.57 


27.82 


- 


25.12 


.60 


14.155 


- 


9.965 


19.15 


- 


11.25 


25.55 


- 


19.12 


.65 














25.54 


- 


16.05 


.7 


12.051 


- 


6.851 


17.58 


- 


8.124 


21.75 


- 


15.57 


.8 


10.119 


- 


4.890 


15.94 


- 


5.994 


19.25 


- 


9.958 


.9 


8.467 


- 


5.616 


14.21 


- 


4.505 


17.15 


- 


7.4i 9 


1.0 


7.110 


- 


2.765 


12.50 


- 


5.451 


15.19 


- 


5.625 


1.1 


6.018 


- 


2.175 


10.87 


- 


2.695 


15.54 


- 


4.554 


1 .2 


5.148 


- 


1.758 


9.40 


- 


2.15 


11.68 


- 


5.591 


1 .5 


4.455 


- 


1.454 








10.18 


- 


2.695 
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TABLE I (Continued) 



Cu + 






Ag + 



Au + 



r 


z 




~u 


Z 




-v" 


Z 




' ir ~ 


1.4 


5.896 




1 .228 


7.03 


- 


1.47 


8.88 


- 


2.14 


1.5 














7.76 


- 


1.78 


1 .6 


5.070 


- 


.9275 


5.51 


- 


1 .08 


6.78 


- 


1.48 


1.8 


2.507 


- 


.7587 


4.11 


- 


.835 


5.25 


“ 


1 .08 


2.0 


2.107 


- 


.6132 


3.26 


- 


.674 


4.14 




.825 


2.2 


1.817 


- 


.5256 


2.65 


- 


.566 


5.54 


- 


.638 


2.4 


1 .606 


- 


.4617 


2.21 




.490 


2.12 


- 


.547 


2.6 


1 .448 


- 


.4135 


1.89 


- 


.434 


2.28 


- 


.468 


2.8 


1 .550 


- 


.5755 


I.65 


- 


.591 


1.95 


- 


.411 


5.0 


1.24J 


- 


.5451 


1.48 


- 


.558 


1.70 


- 


.368 


5.2 


1.179 


- 


.3201 


1.55 


- 


.551 


1.52 


- 


.555 


5.4 


1.150 


- 


.2990 


1.25 




.505 


1.58 


- 


.509 


5.6 


1.095 


- 


.2808 


1.19 




.285 


1 .28 


- 


.287 


5.8 


1 .068 


- 


.2651 


1.15 


- 


.267 


1.21 


- 


.269 


4.0 


1.050 


- 


.2512 


1.10 




.252 


1.15 


- 


.255 


4.5 


1 .022 


- 


.2226 


1 .04 


- 


.225 


1.07 


- 


.223 


5.0 


1 .008 


- 


.2001 


1 .02 


- 


.200 


1.03 


- 


.200 


5.5 


1.003 


- 


.1818 


1 .01 


- 


.182 


1 .01 


- 


.182 


6.0 


1 .001 


— 


.1667 


1.00 


— 


.167 


1 .00 




.167 
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go that (l) becomes 



(19) 



d ^ 



h^( 6 w)e= o 

-n 7 - 



If E is in Rydberg units we can write 



( 20 ) 



■E ^ 



r^i <2. 4 , 



5 



where £ is a dimensionless number, since a Rydberg unit is 
given by 



( 21 ) 



Also, if 



Eclo - 



- A 
r~~\ d 

T^* 



( 22 ) 



or r 



P 






then 



(23) 



v -- 




Inserting (20) and (23) in (19), 



(24) 



But, 









(25) 



so (24) becomes 






AL 

(W) 



(26) 



^ 1- I £ - ■e(f)]R=0 
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where 



(27) 



<8=2. 



or 



Equation (26) can be solved readily using a Runge-Kutta 
4th order approximation. Let 

< 28 ) ^ = $$ 



Then 

(29) 






' d.p’- 

and the set of equations to be solved is 

<» il-t-if'-.i, 

Vdp ' ^ 

The formulas for the solution are 

Z A ,n/ r) 

(51) | A J 1 



where the quantities on the right are calculated in the following 
order: 

^ ^ ^ ,.^4- £ ^ 

A'"(^~ 4'Cfo' 1 ' t '^) K 

= j Q>U, )k 



(52) 



Here, 



(55) 



* 

g-= Wtf)-€.]R 

^ f,-Po 
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After making the calculations 02) and inserting in 
formulas 01 ) we get 

j-i <e(ktk) e,K 



■)-e] 



( 5 4) * r - 4 ^(f 8 )g»-ue t+ -4^(£2i&)e,+^(^) A .K 

- 3 & ^ *- £iK : ) - 4 ] [d?(f 0 ) - e] 4- do ( <f , > e „ 

* * v h>(0-fel[«(.^ -ele. + i 

As before, \ ' — f - ) is calculated using a Gregory- 

Newton second degree divided difference interpolation formula 
of the same form as equation (l6). 

Since equations (jk) give only successive increments in 
R and y it is necessary to obtain starting values for each. 
Assume that near the origin 

(55) !oo-k , 

Also, for small £ , 

( 56 ) 4= - A-f- 

V 

v;here Z is the complete charge on the ion c Differentiating 
(55) twice and inserting it and ( 56 ) in ( 26 ) gives 

^ <*c> 1 2-^lr> { ) + C t> b^t- € b>| 4-2- 3: bi.) ^ 

( 57 ) 4 -C>i^+Gt v 4z-x^)e v tGob5-t£b^75rj 54 v 3 



Thus, 

(58) 



4 - ILb. - 0 



4 - 0 

JO 2. = ^ JO, 

b 5 = -fc 

jOX. • 



12 



For small 



, we want 

(59) K v ^ ^ ^ P 

Inserting (59) in (55) 



(4o) R(6^) = | J( (6^') + U(‘n^) V 4- . • 




so that b, - i • 

As starting values, let us calculate initial values for 
R and y at ~p =.002. Inserting the coefficients in (55) 

' z & e- 

(41) 






We can ignore t compared 
values for any £ are 

Cu + 

R .002 

y .887 



with 2Z^. Then the starting 



Ag + 


Au + 


.002 


.002 


.821 


.709 



The CDC~15o4 program for the calculation is in Appendix 2. 
The magnitude of the initial energy was set at .4 in each case 
and increased in increments of .01 to 1.0 (to 1 .6 for gold). 

The potentials for =.002 were calculated from equation (5 8 ). 
Potentials for (°> 6 were calculated from 



(42) 



\ 




The wave functions from this calculation are intended 

merely for starting values in equation (5), so no great accuracy 

8 i 5 

is required. The error is of the order of h . Only a sample 



^ICunz, K., op. 



cit . , 



pg 188 
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curve of R vs is shown for each metal, that for the energy 
value at which the minimum of the curve oft vs was achieved. 
This is given in Graph 2, The curves of c vs are shown in 
Graph J» 

To solve (j) it may first be put in dimensionless form: 



'-O 



(4?) 



where 



(44) 



and 



w 



* [ e . d ^ ^ (• (f , r) C (r) \< -o 



'.(frbl 2 4*,, c (pi 4>,j (>4 









. <-1 
4-1 



■Jl 



r p 

'~>P 



Here the f . refer to the wave functions for the individual 

n, 1 

electronic shells. These are available (cf note 4). 

For Cu + , the filled shells are the Is, 2s, 2p, Js, 5?> and 
Jd, corresponding to n,l values of 10, 20, 21, JO, J1 , J2. The 
integral then becomes (using r ^ as an example) 

I Kp.rUGMr r 2^? C(r1 1 r + 

{h6) + U. yp V ^ ,(rH?(rUf 

*• ^ I/ 5 ,,(.dC(rUr + jhS J/hj ^ ^ 

and similarly for r 



All the integrals on the right can be evaluated using 
numerical methods with the wave functions cal< ulated above as 
starting values, although the programming will be tedious. 
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The integration, for all practical purposes, probably need not 
proceed beyond the range 0 to 6 , The right hand side of ( 46 ) 



then becomes 




where the a's are the constants resulting from the integrations. 
It is then necessary to solve the set of equations 



by a Runge-Kutta approximation using equations ) and (^2) 
as before. The resulting wave functions are then used to 
re-evaluate the integrals on the right of (45) and so on. The 
procedure is terminated when succeeding sets of wave functions 
differ by less than the required degree of accuracy. 

Unfortunately, there was not enough time available to write 
the CDC-16o4 program. 
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CONCLUSIONS 



The results obtained differ significantly fron thoseof other 
investigators. For example, Fuchs calculated the minimum energy for 
copper as approximately *8, occurring at about 2.9 Bohr radii. 

He gives no estimate of his possible error, however. 

Since the error in finding each increment in the method 
used here is of the order of h , the cumulative error in R to the 
radius at which the minimum occurs is in no case greater than . 005 « 
The accuracy of the method cannot be called into question, therefore. 

There is always the possibility of a mistake in algebra or 
programming, such as the omission of a factor of two. Every step 
in the calculation has been checked over four tines, and if there is 
such an error it must be very subtle in nature. 

It should not be thought, however, that the results of 
previous investigators in this field are necessarily correct. They 
were calculating by hand, and the methods which they were forced to 
use were so oversimplified that reasonable accuracy was often 
impossible. Indeed, if the previous investigators (who all seem 
to have used the same methods) had been in error only .01 units 
for each increment of R, their total error could be very large. 
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APPENDICES 



ITote: figures in the 2 Ci column refer to Aj_ + 
and in the column to Au + , where different 
from the 1 Su column, which is for Cv + 
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APPENDIX I 





ORG 




10000 


STOR 


EQU 




11000 


ANS 


EQU 




12000 


VR1 


EQU 




11100 




LDA 


0 


VIN 




STA 


0 


VR1 


SFT 


ENI 


1 


00000 




ENI 


2 


00000 


BEGIN 


LDA 


0 


ALPHA 




FI-1U 


1 


RAD+1 




STA 


0 


STOR 




LDA 


0 


BETA 




FLRJ 


1 


RAD+2 




FAD 


0 


STOR 




STA 


0 


STOR 




FDV 


0 


TWO 




STA 


0 


STOR+1 




LDA 


0 


STOR 




FMU 


0 


STOR 




STA 


0 


STOR 




LDA 


0 


BETA 




FIIU 


1 


RAD+1 




STA 


0 


STOR+2 




LDA 


0 


ALPHA 




FKU 


1 


RAD+2 
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FAD 


0 


STOR+2 


STA 


0 


ST OR +2 


FDV 


0 


TWO 


STA 


0 


3TOF+5 


LDA 


0 


STOR+2 


FMU 


0 


STOR+2 


STA 


0 


STOR+2 


LDA 


1 


RAD+1 


FAD 


1 


RAD+2 


STA 


0 


STOR+4 


FDV 


0 


TWO 


STA 


0 


STOR+5 


LDA 


0 


STOR+4 


FMU 


0 


STOR+4 


STA 


0 


stor+4 


LDA 


1 


ZSTA 


FSB 


1 


ZSTA+1 


STA 


0 


STOR+5 


LDA 


1 


ZETA+1 


FSB 


1 


ZETA+2 


STA 


0 


STOR+7 


LDA 


1 


RAD 


FSB 


1 


RAD+1 


STA 


0 


STOR+8D 


LDA 


1 


RAD+1 


FSB 


1 


RAD+2 


STA 


0 


STOR+9D 
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LDA 


1 


RAD 


FSB 


1 


RAD+2 


STA 


0 


STOR+1 OD 


LDA 


0 


3T0R+6 


FKU 


0 


ST0R+9D 


STA 


0 


ST0R+11D 


LDA 


0 


STOR+7 


FKU 


0 


ST0R+8D 


STA 


0 


3T0R+12D 


LDA 


0 


STOR+1 ID 


FSB 


0 


3T0R+12D 


STA 


0 


3T0R+1 ID 


LDA 


0 


STOR+CD 


FKU 


0 


3TOR+9D 


FKU 


0 


3T0R+1 OD 


STA 


0 


3T0R+9D 


LDA 


0 


3TOR+1 


FSB 


1 


RAD 


STA 


0 


STOR+1 JD 


FBV 


0 


ST0R+8D 


FKU 


0 


STOR+o 


STA 


0 


STOR+1 4D 


LDA 


0 


STOR+1 


i Ojj 


1 


P* D+1 


FIIU 


0 


STOR+1 pD 


FDV 


0 


STOR+1 OD 


FKU 


0 


STOR+1 ID 


FAD 


0 


STOR+1 4D 


FAD 


1 


ZETA 


FDV 


0 


STOR 


STA 


0 


STOR+1 



22 



LDA 


0 


STOR+5 


FSB 


1 


RAD 


STA 


0 


STOR+1 5D 


FDV 


0 


ST0R+8D 


FMU 


0 


STOR+5 


STA 


0 


STOR+1 6d 


LDA 


0 


STOR+5 


FSB 


1 


RAD+1 


FMU 


0 


STOR+1 5D 


FDV 


0 


STOR+1 OD 


FMU 


0 


STOR+1 ID 


FAD 


0 


STOR+1 6 D 


FAD 


1 


ZETA 


FDV 


0 


ST OR +2 


FAD 


0 


STOR+1 


FLU 


0 


FIVE 


STA 


0 


STOR+5 


STA 


0 


STOR+5 


LDA 


0 


STOR+5 


FSB 


1 


RAD 


STA 


0 


STOR+1 7D 


FDV 


0 


3T0R+8D 


FMU 


0 


stor+6 


STA 


0 


STOR+1 8D 


LDA 


0 


STOR+5 


FSB 


0 


RAD+1 



CQNV 



FMU 0 STOR+17D 
FDV 0 STOR+IOD 
FMU 0 3T0R+11D 
FAD 0 ST0R+18D 
FAD 1 ZSTA 
FDV 0 STOR+4 
FMU 0 FIGHT 
FAD 0 3T0R+5 
FMU 0 TWO 
FDV 0 NIKE 
ST A 0 3T0R+3 
LDA 1 RAD+2 
FSB 1 RAD+1 
FMU 0 ST OR +5 
FAD 0 VR1 
STA 0 VR1 
LDQ 1 RAD+2 
3TQ 2 AITS 
STA 2 AN3+1 
INI 2 00002 
ENI 0 00000 

ISK 1 55D ISK 1 59D ISK 1 49D 

5LJ 0 BEGIN 

ENI 5 00000 

ENI 4 ANS 

slj 4 71000 
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E.TI 


0 00000 


001 


0 ANS 


006 


0 0 


INI 


A 00006 


SNA 


A 00000 


SAU 


0 CONV+1 


ENI 


0 00000 


ISK 


5 17D ISK 5 19D ISK 5 16 D 


SLJ 


0 CONV 


ENI 


0 00000 


SL3 


0 00000 


VIN DSC 


-.18181818181818 


ALPHA DEC 


1 . 77 A 596666 


BETA DEC 


.225A055A 


TWO DEC 


2. 


FIVE DSC 


5 . 


EIGHT DEC 


8. 


NINE DSC 


9. 


RAD DEC 


6. 


DSC 


5.5 


DEC 


5 . 


DEC 


4.5 


DEC 


A. 


DEC 


5.8 


DEC 


5.8 


DEC 


5. A 
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DSC 


5.2 










DSC 


5. 










DEC 


2.8 










DSC 


2.6 










DEC 


2.4 










DEC 


2.2 










DSC 


2. 










DEC 


1.8 










DEC 


1.6 










DEC 


1.4 


DSC 


1.4 


DSC 


1.5 


DSC 


1.5 


DSC 


1.2 


DEC 


1 .4 


DEC 


1.2 


DSC 


1.1 


DSC 


1.5 


DEC 


1.1 


DSC 


1.0 


DSC 


1.2 


DSC 


1. 


DSC 


.9 


DEC 


1.1 


DSC 


.9 


DEC 


.8 


DEC 


1. 


DEC 


.8 


DEC 


.7 


DEC 


.9 


DEC 


.7 


DEC 


.6 


DEC 


.8 


DEC 


.6 


DEC 


.55 


DSC 


.7 


DEC 


.55 


DEC 


.5 


DSC 


.85 


DEC 


.5 


DEC 


.45 


DSC 


.6 


DSC 


.45 


DEC 


.4 


DSC 


.55 


DEC 


.4 


DEC 


.58 


DEC 


.5 


DSC 


• 55 


DEC 


• 56 


DEC 


.45 


DEC 


• 5 


DEC 


.54 


DEC 


.4 


DSC 


.28 


DSC 


• 52 


DSC 


.58 


DSC 


.26 


DSC 


.5 


DEC 


.58 



2 6 



ZETA 



DEC 


.24 


DEC 


.28 


DSC 


• 5 Zf 


DEC 


.22 


DSC 


.26 


DEC 


.52 


DEC 


.2 


DEC 


.24 


DEC 


.5 


DEC 


.18 


DEC 


.22 


DEC 


.28 


DSC 


.16 


DEC 


.2 


DEC 


.26 


DEC 


.14 


DEC 


.18 


DEC 


.24 


DEC 


.12 


DEC 


.16 


DEC 


.22 


DEC 


.1 


DEC 


.14 


DEC 


.2 


DEC 


.09 


DEC 


.12 


DEC 


.18 


DEC 


.08 


DEC 


.11 


DSC 


.16 


DEC 


.07 


DEC 


.1 


DSC 


.14 


DEC 


.06 


DEC 


.09 


DEC 


.12 


DEC 


.05 


DEC 


.08 


DEC 


.1 


DEC 


.04 


DEC 


.07 


DEC 


.08 


DEC 


.055 


DEC 


.06 


DEC 


.06 


DEC 


.05 


DEC 


.05 


DEC 


.04 


DEC 


.025 


DEC 


.04 


DEC 


.02 


DEC 


.02 


DEC 


.055 


DEC 


0 


DEC 


.015 


DSC 


.05 


ZETA DSC 


1. 


DEC 


.01 


DSC 


.025 


DEC 


1.01 


DEC 


.005 


DSC 


.02 


DEC 


1.05 


DEC 


0 


DEC 


.015 


DSC 


1.07 


DSC 


1 .001 


DSC 


.01 


DEC 


1.15 


DEC 


1.005 


DEC 


.008 


DEC 


1.21 


DEC 


1.008 


DSC 


• 

O 

O 

On 


DEC 


1.28 


DSC 


1 .022 


DSC 


.004 


DEC 


1.58 
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DSC 


1.050 


DEC 


.002 


DEC 


1.52 


DEC 


1.068 


DSC 


0 


DEC 


1.70 


DEC 


1.095 


ZETA DEC 


1 


DEC 


1.95 


DEC 


1.150 


DEC 


1.01 


DSC 


2.28 


DEC 


1.179 


DEC 


1.02 


DSC 


2.72 


DEC 


1.245 


DEC 


1.04 


DEC 


5.5^ 


DEC 


1.550 


DEC 


1.10 


Di^C 


4.14 


DEC 


1 .448 


DEC 


1.15 


DEC 


5.25 


DEC 


1 .606 


DEC 


1.19 


DEC 


6.78 


DSC 


1.817 


DEC 


1.25 


DEC 


7.76 


DEC 


2.107 


DEC 


1.55 


DSC 


8.88 


DEC 


2.507 


DEC 


1.48 


DEC 


10.18 


DSC 


5.070 


DEC 


1.65 


DEC 


1 1 .68 


DEC 


5.896 


DSC 


1 .89 


DEC 


15.54 


DEC 


4.455 


DEC 


2.21 


DEC 


15.19 


DSC 


5.148 


DEC 


2.65 


DEC 


17.15 


DEC 


6.018 


DEC 


5.26 


DEC 


19.25 


DEC 


7.110 


DEC 


4.11 


DSC 


21.75 


DEC 


8.457 


DEC 


5.51 


DEC 


25.54 


DSC 


10.119 


DEC 


7.05 


DEC 


25.55 


DSC 


12.051 


DEC 


9.4 0 


DEC 


27.82 


DSC 


14.155 


DEC 


10.87 


DEC 


50.96 


DEC 


15.196 


DEC 


12.5 


DEC 


54.75 


DEC 


16.195 


DEC 


14.21 


DEC 


58.94 


DEC 


17.114 


DSC 


15.94 


DEC 


40.65 


DEC 


17.964 


DEC 


17.58 


DEC 


42.56 
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DSC 


18.791 


DEC 


19.15 


DEC 


44.05 


DEC 


19.756 


DEC 


20. 04 


DEC 


45.62 


DEC 


20.187 


DEC 


21.16 


DEC 


47.15 


DEC 


20.699 


DEC 


22.65 


DEC 


48.62 


DEC 


21.279 


DEC 


24. 61 


DEC 


50.09 


DEC 


21 .929 


DEC 


25.55 


DEC 


51.66 


DSC 


22 .646 


DEC 


26.57 


DEC 


55.40 


DEC 


25.415 


DEC 


27.66 


DEC 


55.44 


DEC 


24.209 


DEC 


28.8 


DEC 


57.75 


DEC 


24.977 


DEC 


29.98 


DSC 


60.51 


DSC 


25.689 


DEC 


51.16 


DEC 


62.92 


DEC 


26.520 


DEC 


52.51 


DEC 


65.41 


DEC 


26.606 


DEC 


55*45 


DEC 


67.76 


TM ?n 
i 


26.884 


DEC 


54.48 


DEC 


70.14 


DEC 


27.161 


DEC 


55.^8 


DEC 


72.86 


DEC 


27.451 


DSC 


56.48 


DEC 


75.55 


DSC 


27.765 


DEC 


57.55 


DSC 


77.55 


DEC 


28.106 


DEC 


58.71 


DEC 


79. 


DEC 


28.280 


DEC 


4o.o4 


END 




DEC 


28.451 


DEC 


40.78 






DSC 


28.615 


DEC 


41.52 






DEC 


28.759 


DEC 


42.26 






DEC 


28.875 


DSC 


42.96 






DEC 


26.948 


DSC 


45.62 






DEC 


28.992 


DEC 


44.20 






DEC 


29.000 


DSC 


44.75 






END 




DEC 


45.25 







29 



DEC 



45-5 

DEC 45.78 
DEC 46.06 
DEC 46.56 
DEC 46.65 
DEC 46.85 
DEC 46.91 
DEC 46 . 95 
DEC 46.99 
DEC 47 . 

DSC 47 . 

END 

Note 1 : Library subroutine DECOF 1 must be placed in memory 

starting at cell 71000. Output will be on magnetic 
tape unit 4. 

Note 2: Program prints 0 and 'V . Two to a line in decimal 

floating point. 
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APPENDIX II 





REM 




SCKR EQ BY RK METHOD 




ORG 




20000 


7ALUE 


EQU 




22000 


3T0R 


EQU 




21000 


ViYE 


EQU 




21600 


ARE 


EQU 




21650 


EilERC-Y 


EQU 




21700 




LDA 


0 


IITEN 




STA 


0 


ENERGY 




Eli I 


3 


00000 




EITI 


0 


00000 


BSC- II! 


El! I 


i 


00000 




Eli I 


2 


00000 




LDA 


0 


rut 




STA 


0 


ARE 




LDA 


0 


YIN 




STA 


0 


WYE 


RETURN 


LDA 


1 


RAD 




FSB 


1 


RAD+1 




STA 


0 


3T0R 




LDA 


1 


RAD 




FSB 


1 


RAD+2 




STA 


0 


STOR+1 




LDA 


1 


RADhI 




FAD 


1 


RADh 2 
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•- 



FSB 


1 


RAD 


FSB 


1 


RAD 


STA 


0 


STOR+2 


LDA 


1 


POT+1 


FSB 


1 


POT +2 


FMU 


0 


STOR+2 


FDV 


0 


FOUR 


FDV 


0 


STOR+1 


STA 


0 STOR+5 


LDA 


1 


POT 


FSB 


1 


POT+1 


FMU 


0 


STOR+2 


FMU 


0 


STOR+2 


FDV 


0 


FOUR 


FDV 


0 


STOR 


FDV 


0 


STOR+1 


STA 


0 


STOR+4 


LDA 


1 


POT 


FSB 


0 


STOR+4 


FAD 


0 STOR+5 


STA 


0 


STOR+5 


LDA 


1 


RAD+1 


FSB 


1 


RAD 


STA 


0 STOR+6 


FMU 


0 


3T0R+5 


STA 


0 


ST OR +7 
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FMU 


0 


stor+6 


3TA 


0 


ST0R+8D 


LDA 


1 


POT 


FSE 


0 


ENERGY 


STA 


0 


STOR+9D 


LDA 


0 


STOR+5 


FSB 


0 


ENERGY 


FMU 


0 


STOR+9D 


FiU 


0 


3T0R+8D 


flu 


0 


ARE 


FDV 


0 


FOUR 


STA 


0 


3T0R+10D 


LDA 


0 


ENERGY 


FMU 


0 


WYE 


FMU 


0 


3TOR+7 


STA 


0 


STOR+1 ID 


LDA 


0 


3TOR+5 


FIU 


0 


WYE 


FMU 


0 


STOR+7 


STA 


0 


STOR+1 2D 


LDA 


0 


3TOR+5 


FIU 


0 


ARE 


FMU 


0 


STOR+6 


FMU 


0 


TWO 


STA 


0 


STOR+1 JD 


LDA 


0 


ENERGY 



FKU 


0 


THREE 


FMU 


0 


ARE 


fi:u 


0 


STOR+o 


STA 


0 


ST0R+14D 


LDA 


1 


POT 


FMU 


0 


ARE 


FMU 


0 


STOR+6 


STA 


0 


ST0R+15D 


LDA 


0 


V/YE 


FI-.U 


0 


SIX 


FAD 


0 


ST0R+15D 


FSB 


0 


ST0R+14D 


FAD 


0 


ST0R+15D 


FAD 


0 


5T0R+ 1 2D 


F3E 


0 


ST0R+11D 


FAD 


0 


ST0R+10D 


FMU 


0 


STOR+6 


FDV 


0 


SIX 


STA 


0 


STOR+16D 


LDA 


1 


POT+1 


FMU 


0 


LYE 


FMU 


0 


STOR+o 


STA 


0 


STOR+1 8D 


LDA 


0 


VYE 


FI-IU 


0 


stor+5 


FDV 


0 


TV, r O 
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PAD 


0 


ARE 


STA 


0 


3TCR+19D 


LDA 


0 


STOR+5 


FSB 


0 


ENERGY 


STA 


0 


ST0R+20D 


LDA 


1 


POT+1 


FSB 


n 


ENERGY 


FMU 


0 


ST0R+20D 


FMU 


0 


STOR+1 9D 


FMU 


0 


3T0R+7 


FDV 


0 


TWO 


STA 


0 


STOR+21 D 


LDA 


0 


ENERGY 


FMU 


0 


ARE 


FI IU 


0 


SIX 


STA 


0 


ST0R+22D 


LDA 


1 


POT+1 


FLU 


0 


ARE 


STA 


0 


STOR+25D 


LDA 


1 


POT 


FSB 


0 


ENERGY 


STA 


0 


STOR+24D 


LDA 


0 


STOR+5 


FSB 


0 


ENERGY 


fi:u 


0 


STCR+24D 


FMU 


0 


ARE 
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FMU 


0 


STOR+7 


FDV 


0 


TUC 


3TA 


0 


ST0R+25D 


LDA 


0 


ETERC-Y 


FI-X- 


0 


’.TE 


?i:u 


0 


sroR+6 


FMU 


0 


THREE 


3TA 


0 


ST0R426D 


LDA 


0 


STOR+5 


FMU 


0 


V/YE 


FMU 


0 


3T0R+5 


FMU 


0 


tv:o 


r>rr\ f 
Oi i 1 


0 


ST OR +27 D 


LDA 


0 


3T0R+5 


FMU 


0 


ARE 


FMU 


0 


FOUR 


3TA 


0 


ST0R+2GD 


LDA 


1 


POT 


Fx-IU 


0 


ARE 


F3B 


0 


3T0R+22D 


r ID 


0 


3T0R+28D 


FAD 


0 


3T0R+27D 


FSB 


0 


ST0R+26D 


FAD 


0 


3T0R+25D 


FAD 


0 


S"OR+25D 


FAD 


0 


3T0R+21D 


FAD 


0 


3TOR+18D 



CONV 



PI-IU 0 STCR+6 
FDV 0 SIX 
FAD 0 WYE 
STA 0 WYE 
3TA 2 VALUE+5 
LDA 1 RAD+1 
STA 2 VALUE 
LDA 0 STOR+1 5D 
FAD 0 ARE 
STA 0 ARE 
STA 2 VALUE+1 
FDV 1 RAD+1 
STA 2 VALUE+2 
INI 2 00004 
ENI 0 00000 

ISK 1 56D I5K 1 59U ISK 1 52D 

SLJ 0 RETURN 

ENI 4 00000 

ENI 5 VALUE 

SLJ 4 71000 

ENI 0 COOOO 

001 0 VALUE 

004 0 0 

INI 5 00004 

ENA 5 00000 

3AU 0 COXV+1 
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EKI 


0 00000 












ISK 


4 56 D 


ISK 


4 59D 


ISK 


4 52D 




3LJ 


0 COITV 












LDA 


0 ENERGY 












FAD 


0 DELES 












STA 


0 ENERGY 












EKI 


0 ooooc 












ISK 


5 00075 






ISK 


5 00170 




SLJ 


0 BEGIN 












EKI 


0 00000 












SL3 


0 00000 










IKSN 


DSC 


-.4 










TV.’O 


DEC 


2. 










THREE 


DEC 


5. 










FOUR 


DEC 


4. 










SIX 


DEC 


6. 










DELEFT 


DEC 


-.01 










RIK 


DSC 


.002 










YE 


DEC 


.887 


DEC 


.821 


DEC 


.709 


RAD 


DSC 


.002 












DEC 


ITS 

O 

O 

« 


DEC 


-3- 

o 

o 

• 


DEC 


.02 




DEC 


.01 


DSC 


.006 


DEC 


.04 




DEC 


.015 


DEC 


.008 


DEC 


.06 




DEC 


.02 


DEC 


.01 


DEC 


.08 




DSC 


.025 


DEC 


.015 


DEC 


.1 




DEC 


to* 

o 

« 


DEC 


.02 


DEC 


.12 
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DEC 


.055 


DEC 


.025 


DEC 


.14 


DEC 


.04 


DEC 


.05 


DEC 


. 16 


DEC 


.05 


DEC 


.055 


DEC 


.18 


DEC 


.06 


DEC 


.04 


DEC 


.2 


DEC 


.01 


DEC 


.05 


DEC 


.22 


DEC 


.08 


DEC 


.06 


DEC 


.24 


DSC 


.09 


DEC 


.07 


DEC 


.26 


DEC 


.1 


DSC 


.08 


DEC 


.28 


DSC 


.12 


DEC 


.09 


DSC 


.5 


DEC 


.14 


DEC 


.1 


DEC 


.52 


DEC 


.16 


* DEC 


.11 


DEC 


• 54 


DEC 


.18 


DEC 


.2 


DSC 


.58 


DSC 


.2 


DEC 


.14 


DSC 


.58 


DEC 


.22 


DSC 


.16 


DEC 


.4 


DEC 


.24 


DSC 


.18 


DSC 


.45 


DEC 


.26 


DSC 


.2 


DEC 


.5 


DEC 


.28 


DEC 


.22 


DEC 


.55 


DEC 


.5 


DEC 


.24 


DEC 


.6 


DEC 


.55 


DEC 


.26 


DEC 


.65 


DEC 


.4 


DEC 


.28 


DEC 


.7 


DSC 


.45 


DEC 


.5 


DEC 


.8 


DEC 


.5 


DSC 


.52 


DEC 


.9 


DEC 


.55 


DEC 


.54 


DEC 


1 . 


DEC 


.6 


DSC 


.56 


DEC 


1.1 


DSC 


.7 


DEC 


.58 


DEC 


1.2 


DEC 


.8 


DEC 


.4 


DEC 


1.5 
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DEC 


.9 


DSC 


.45 


DEC 


1.4 


DEC 


1 . 


DSC 


.5 


DSC 


1.5 


DEC 


1.1 


DSC 


.55 


DSC 


1.6 


DEC 


1.2 


DSC 


.6 


DSC 


1.8 


DEC 


1.5 


DSC 


.7 


DSC 


2. 


DSC 


1.4 


DSC 


.8 


DEC 


2.2 


DEC 


1 .6 


DEC 


.9 


DSC 


2.4 


DEC 


1.8 


DEC 


1. 


DSC 


2.6 


DSC 


2. 


DEC 


1.1 


DSC 


2.8 


DSC 


2.2 


DEC 


1.2 


DEC 


5. 


DEC 


2.4 


DSC 


1 .4 


DSC 


5.2 


DEC 


2.6 


DSC 


1 .6 


DEC 


5.4 


DSC 


2.8 


DSC 


1 .8 


DEC 


5.6 


DEC 


5. 


DEC 


2. 


DSC 


5.8 


DEC 


5.2 


DSC 


2.2 


DEC 


4. 


DSC 


5.4 


DEC 


2.4 


DEC 


4.5 


DEC 


5-6 


DEC 


2.6 


DSC 


5. 


DSC 


5.8 


DEC 


2.8 


DSC 


5.5 


DEC 


4. 


DSC 


5. 


DSC 


6 


DEC 


4.5 


DEC 


5-2 


DEC 


6.5 


DEC 


5. 


DSC 


5.4 


DEC 


7. 


DEC 


5.5 


DSC 


5.6 


POT DSC 


-79000 


DEC 


6 . 


DSC 


5.8 


DSC 


- 7004 


DSC 


6.5 


DEC 


9 


DSC 


- 51 46, 


DEC 


7. 


DEC 


4.5 


DEC 


- 1907 


DEC 


-29000. 


DEC 


5- 


DEC 


- 1511 



4o 



DEC 


- 11336 . 


DEC 


5.5 


DSC -966.2 


DSC 


- 5539.4 


DEC 


6 . 


DEC -743=8 


DEC 


- 5610 . 


DEC 


6.5 


DEC -590.6 


DEC 


- 2648.8 


DSC 


7. 


DEC -480.6 


DEC 


- 2074.6 


POT DEC 


- 27820 . 


DEC -398.8 


DSC 


- 1694.1 


DEC 


-14880. 


DSC -336. 


DEC 


- 1423.9 


DEC 


- 9974. 


DSC -286.6 


DSC 


- 1222.5 


DSC 


7410. 


DEC -246.8 


DEC 


- 943.2 


DSC 


- 5858. 


DSC -214.2 


DSC 


- 759.16 


DSC 


- 3890. 


DEC -187.1 


DEC 


- 629.14 


DEC 


- 2876 . 


DEC -164,5 


DEC 


- 552.6 


DEC 


- 2258 . 


DEC -144.9 


DEC 


- 458.28 


DEC 


- 1 S4i . 


DSC -128.4 


DSC 


- 399.42 


DEC 


- 1545. 


DSC -114.3 


DEC 


- 512.54 


DSC 


- 1320 . 


DEC -102.2 


DSC 


- 252.1 


DSC 


- 1019. 


DEC - 91.68 


DEC 


- 208 . 14 


DSC 


- 818.4 


DSC - 71.3 


DSC 


- 175.07 


DSC 


- 675.6 


DSC - 56=82 


DSC 


- 149.49 


DSC 


- 569.4 


DEC - 46.24 


DEC 


- 129.25 


DSC 


- 487.6 


DEC - 38.24 


DSC 


- 112.9 


DEC 


- 423.2 


DEC - 52.06 


DEC 


99.466 


DSC 


- 371.4 


DEC - 27.14 


DSC 


88.248 


DSC 


- 329. 


DEC - 19.88 


DSC 


- 78.758 


DEC 


- 265.6 


DSC - 14.84 


DEC 


50.448 


DEC 


- 219 . 


DEC - 11.25 


DEC 


- 47.5 


DEC 


185.6 


DEC - 8.668 



4i 



DEG 


« 


57.526 


DEC 


- 


155. 


DEC 




6.782 


DEC 


- 


50.100 


DSC 


- 


155.9 


DSC 


- 


5.59 


DEG 




24.578 


DEC 


- 


116.1 


DSC 


- 


4.56 


DSC 


- 


19.926 


DSC 


- 


101 .6 


DEC 


- 


5.56 


DEC 


- 


15.702 


DSC 


- 


89.52 


DEC 


~ 


2.96 


DSC 


- 


9.78 


DEC 


- 


79.5 


DEC 


- 


2.16 


DSC 




7.252 


DSC 




71 .06 


DEC 


- 


I.65 


DSC 


- 


5.526 


DEC 


- 


65.9 


DSC 


- 


1.52 






4.55 


DEC 


- 


57.78 


DSC 


* 


1.09 


DSC 


- 


5.516 


DSC 


- 


52.48 


DSC 


- 


.956 


DSC 


- 


2.908 


DEC 




47.92 


DSC 


- 


.822 


DSC 


“ 


2.456 


DSC 


- 


58.96 


DSC 




.756 


DEC 


- 


1.855 


DSC 




52.12 


DEC 


— 


.6 7 


DSC 


- 


1.477 


DSC 


- 


26.74 


DEC 




.618 


DSC 


- 


1 .226 


DSC 


- 


22.46 


DEC 


~ 


.574 


DSC 


- 


1 .051 


DSC 


- 


16.25 


DEC 


- 


.558 


DSC 


- 


.9254 


DSC 


- 


11.99 


DEC 


- 


.506 


DSC 




.8266 


DSC 


- 


9.01 


DSC 


- 


.446 


DSC 


- 


.7510 


DEC 


- 


6.902 


DEC 


* 


.4 


DEC 


- 


.6902 


DSC 


- 


5.59 


DEC 


- 


.564 


DEC 


- 


.5402 


DEC 


- 


4.5 


DSC 




• 554 


DSC 


- 


.598 


DSC 


- 


2.94 


DSC 


- 


.508 


DEC 


- 


.5616 


DSC 


- 


2.16 


DSC 




.286 


DSC 




.5502 


DSC 


- 


1.666 


END 






DSC 


- 


.5024 


DSC 




1.548 








DSC 


- 


.4452 


DEC 


— 


1.152 









42 



DEC 


- . 4002 


DSC 


- 


.98 


DSC 


- .3636 


DEC 


- 


.868 


DEC 


- . 333 ^ 


DEC 


- 


.782 


DEC 


- .5077 


DSC 


- 


.716 


DEC 


- .2357 


DEC 


- 


.662 


END 




DEC 


- 


.61 






DEC 


* 


.5 7 






DEC 


- 


.554 






DEC 


- 


.504 






DEC 


- 


.446 






DEC 


- 


.4 






DEC 


- 


.384 






DEC 


- 


.554 






DSC 


- 


.308 






DEC 


- 


.286 






END 







Note 1 : Note 1 to Appendix I also applies here 

Note 2: Program prints p t? ^ ^ in decimal floating point 

’ e d.(> 

along each line. 
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